Conical intersections in laboratory coordinates with ultracold molecules 
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For two states of opposite parity that cross as a function of an external magnetic field, the addition of an 
electric field will break the symmetry and induce an avoided crossing. A suitable arrangement of fields may be 
used to create a conical intersection as a function of external spatial coordinates. We consider the effect of the 
resulting geometric phase for ultracold polar molecules. For a Bose-Einstein condensate in the mean-field ap- 
proximation, the geometric phase effect induces stable states of persistent superfluid flow that are characterized 
by half-integer quantized angular momentum. 
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It is well known that the potential energy surfaces for 
molecular electronic states of the same symmetry can cross 
at a point in 2 dimensions or on a surface of dimension n — 1 
in n dimensions. These crossings are known as conical in- 
tersections because the two surfaces locally form a double 
cone. Conical intersections have a wealth of interesting conse- 
quences for molecular structure and dynamics. For example, 
they are responsible for the Jahn-Teller effect UJJ and play an 
important role in nonadiabatic processes [2]. One of the most 
interesting consequences of conical intersections is the geo- 
metric phase (Berry phase) effect y|]: when the nuclei follow 
a path that encircles a conical intersection once and returns 
to the original configuration (pseudorotation), the electronic 
wavefunction changes sign. Since the total wavefunction must 
be a single-valued function of coordinates, this requires that 
the wavefunction for nuclear motion must also change sign. 
This has important dynamical consequences: it produces half- 
odd-integer quantization for free pseudorotation Hi and may 
have significant effects on collision cross sections ]5[]. 

The purpose of the present paper is to explore conical in- 
tersections of a different type. It is now possible to produce 
atomic and molecular Bose-Einstein condensates (BECs) and 
to subject them to applied magnetic and electric fields. The 
atomic and molecular states split and shift as a function of the 
magnetic field (Zeeman effect) and electric field (Stark effect). 
In the absence of an electric field, parity is conserved, so it is 
possible to tune the magnetic field so that two levels of differ- 
ent parity are exactly degenerate with one another. However, 
if a simultaneous electric field is applied, the two levels of 
opposite parity are mixed and the degeneracy is resolved [6]. 
Conical intersections can thus occur at points where the elec- 
tric field is zero. It is possible to envisage an arrangement of 
fields that creates conical intersections between two atomic or 
molecular levels as a function of external spatial coordinates 
rather than internal coordinates. 

All the results that apply to conical intersections between 
molecular potential energy surfaces will continue to apply in 
this new situation. The internal (electronic/vibrational/spin) 
wavefunction of the atom or molecule will change sign along 
a path that encircles the intersection once, so the spatial wave- 
function of the condensate must also change sign. 

Quantized vortices are a characteristic sign of superfluid 
flow 0], and since their first observation in the context of 



Bose-Einstein condensed dilute atomic gases JpLhave been 
key to some spectacular experimental results |9|,[l0||. Recent 
experiments placing a BEC within a toroidal trapping geome- 
try iTTTt fT^L []~3ll have enabled the observation of persistent flow 
around a toroidal trap 111 111 . The possibility of forming half- 
integer quantized vortices within a spinor atomic BEC config- 
uration has been investigated theoretically lfT^[l5l[l6l[l7ll . as 
have a variety of differing consequences ofgeometric phase 
effects in atomic BEC systems lfl8i [l9l [20T l2lll . In this work 
we combine these different threads to show how, with an ap- 
propriate configuration of magnetic and electric fields, a BEC 
of heteronuclear diatomic molecules will assume a toroidal 
geometry, such that the geometric phase causes the system 
to manifest macroscopically occupied states of half-integer 
quantized persistent flow. 

Effects of this type can in principle be observed in any sys- 
tem where two states of opposite parity can be tuned into de- 
generacy with a magnetic field and can be coupled with an 
electric field. However, for atomic systems states of differ- 
ent parity are usually far apart at zero field. More accessible 
examples are provided by gases of heteronuclear alkali-metal 
dimers such as RbCs and KRb, which are the targets of cur- 
rent experiments. In the present work we illustrate the effect 
for a gas of KRb molecules in a single vibrational level of the 
lowest triplet state, 3 Z + . 

The energy levels of a 3 2 + molecule in an applied field are 
conveniently expanded in a fully decoupled basis set of func- 
tions \NMn)\S Ms ), where N and S are quantum numbers for 
molecular rotation and electron spin and and Ms are the 
corresponding space-fixed projections onto the magnetic field 
axis. Nuclear spin is neglected here for simplicity. A simple 
form of the Hamiltonian that contains all the essential ingre- 
dients is 

H = Bfl 2 +^A(§ 2 -3Sl)+g e fi 3 B z Ms-H.E, (1) 

where B is the molecular rotational constant, A is the spin-spin 
coupling constant, Sz is the projection of S onto the molec- 
ular axis, Bz is the magnetic field orientated along the space- 
fixed Z axis, and E and /j, are the vectors representing the 
electric field and molecular electric dipole moment. KRb has 
not been characterized in detail spectroscopically but ultra- 
cold KRb has recently been formed in the lowest rovibrational 



2 




Magnetic Field Strength (Gauss} 

FIG. 1: (color online). Energy levels of KRb as a function of mag- 
netic field B. The upper panel shows an overview, with each zero- 
field level (labeled by (N,J) splitting into 27+1 components as a 
function of magnetic field. The lower panel shows how two levels of 
different parity cross (Mj = and +1) in the absence of an electric 
field (solid lines) but avoid one another in the presence of a 5 kV/cm 
electric field (dashed lines). 



levels of both singlet and triplet states If22tl . Electronic struc- 
ture calculations give an equilibrium distance r e = 5.901 A 
for the triplet state l23Tl . This allows B and A for the lowest 
vibrational level to be estimated as B — 0.01813 crrT 1 and 
A = -0.00632 cm -1 . The dipole moment function has been 
calculated by Kotochigova, Julienne and Tiesinga [24] and has 
a value around 0.051 D near r e . 

Fig.[T]shows the lowest rotational levels of KRb (a 3 2 + ) as 
function of magnetic field in the presence and absence of a 
moderate electric field (5 kV/cm). At zero field, the N - 
level has a single sublevel with total angular momentum 7 = 
1, while the N = 1 level is split into 3 sublevels with 7 = 0, 
1 and 2. When a magnetic field is applied, each sublevel is 
split into 27+1 components labeled by Mj. At zero electric 
field, the Mj — + 1 level originating from N = 0, 7 = 1 and 
the Mj =0 level originating from N = I, J = I have different 
parity and cross near Bz = 187 G. 

When a non-zero electric field is introduced, parity is no 



longer conserved. However, if the electric and magnetic fields 
are parallel, Mj is conserved and Mj = and +1 states cross. 
We have therefore chosen the electric field to be perpendicular 
to the magnetic field to induce an avoided crossing between 
the Mj — and +1 states as shown in Fig.Q] 

We may envisage an experiment in which a BEC is sub- 
jected to a magnetic field Bz, orientated along the space- 
fixed Z axis, which varies along the X axis with field gradi- 
ent dBz/dX. An inhomogeneous electric field E is oriented 
along the X axis with a magnitude which varies with Y as 
dE/dY, vanishing on a plane at Y = 0. This creates a seam 
of conical intersections along the line 0,0, Z where X = is 
the position at which the magnetic field brings the two states 
into degeneracy. Adding an external cylindrically symmet- 
ric optical trapping potential in the XY plane of the electro- 
magnetic field gradients, V opt (p,(p,Z) = jM(a> 2 p 2 + w|Z 2 ), 
where p 2 — X 2 + Y 2 , creates a toroidally shaped potential 
around the conical intersection, with a radial minimum at po- 
The left-hand side of Figure|2]shows the resulting potential for 
a magnetic field gradient of 5 G/cm, an electric field gradient 
of 6.8 kV/cm 2 and an optical trapping potential with a height 
of 7 pK at p = 30 pm M25I1 centered at the conical intersec- 
tion. This potential has significant anisotropy (about 10 nK), 
which is manifested as an asymmetry along a cut with Y — 
as shown schematically in the lower panel. 

The anisotropy of the toroidal trapping potential V tmp ((f>) 
can be controlled by offsetting the optical trapping poten- 
tial from the point of intersection. Assuming that the KRb 
eigenstates are linear as a function of magnetic field over the 
range of the intersection, with respective gradients a and b, 
the asymmetry along 7 = will be zero when the optical trap 
is centered at xo — (a + b)/(2Ma> p ). The trapping potential 
still has slightly different depths along the X and Y axes, but 
this can be minimized by adjusting the electric field gradient. 
The right-hand side of Figure|2]shows the resulting optimized 
trapping potential, with Xq = 0.0715 pm and dE/dY = 6.723 
kV/cm 2 . This potential has an angular anisotropy on the order 
ofO.OlnK. 

For a molecule of mass M moving in a toroidal potential 
such as those in Fig.|2j the single-particle Schrodinger equa- 
tion is approximately separable, with solutions *P(p, <p, Z) = 
i/r(p)<f>((p)ip(Z). In the absence of anisotropy, the geometric- 
phase-induced anti-periodic boundary condition d>(0) = 
— <&((p + 2n) gives <t>(0) = (27r) -1 ^ 2 exp(z'm0), where m takes 
half-integer values ±A,±|,±|, etc. If the boundary con- 
ditions were periodic over 2n, m would take integer values 
0, +1, +2, etc. For periodic boundary conditions the ground 
state has zero angular momentum m = 0, but for anti-periodic 
boundary conditions it has a nonzero angular momentum m - 
+ \ . The single-particle energy spectrum is 

E = ( VZ + 2~) hC ° Z + ( Vp + 2~) hUp + bmt ~ i) ' (2) 

where the third term represents the rotational energy of the 
particle with a rotational constant b ml = h 2 /2Mp 2 . 

In order for geometric phase effects to be observed, the 
angular anisotropy of the toroidal trapping potential V tmp (4>) 
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FIG. 2: (color online). Toroidal potentials formed around conical intersections. Left-hand panels: potential formed when the optical trapping 
potential (7 nK at 30 /urn) is centered at the point of intersection (p=0), with field gradients 5 G/cm X and 6.8 kV/cm 2 Y. Right-hand panels: 
potential formed when the trapping potential is offset along X by xo = 0.0715 /im, with field gradients 5 G/cm X and 6.723 kV/cm 2 Y. The 
electric field gradient is chosen to minimize the anisotropy in each case. The well depths are given in nK relative to the point of intersection. 
The lower panels show schematic cuts through the potential at Y = 0: KRb eigenstates (red, dashed), optical trapping potential (blue, dot- 
dashed) and the resultant toroidal potential (black, solid). 



must be small enough to allow the wavefunction to fully en- 
circle the intersection. For the potentials shown in Fig. [2] 
jHoj p as 4 nK and b mt ~ 0.5 nK. For the potential on the left- 
hand side, the anisotropy is large compared to b I0t , so that the 
single-particle wavefunction will be localized on one side of 
the trap. However, for the potential on the right-hand side, the 
anisotropy is small compared to b mi and the single-particle 
wavefunction will fully encircle the conical intersection and 
exhibit half-integer quantization. 

A BEC of a dilute gas may be modeled by the Gross- 
Pitaevskii equation (GPE). The time-independent GPE is 

[H + u(<f>)\V(p, 0, Z)| 2 ] ¥(p, 0, Z) = i&(p, 0, Z), (3) 

where is the chemical potential and the mean-field wave- 
function ¥ is normalized to unity, J = 1. The effec- 
tive interaction strength u(4>) is 47rti 2 Na(<f))/M, where N is the 
number of particles in the condensate. The internal molec- 
ular wavefunction at angle <p may be written in terms of the 



individual molecular states \p \ and \p2 as 

Mfi = ifn cos(0/2) + (A 2 sin(0/2), (4) 
so the effective angle-dependent scattering length a(<p) is 



a((p) — | (3fln + 3a22 + 2ai2) + 5(011 - 022) cos < 
+ i (an + a 2 2 - 2ai 2 )cos20, 



(5) 



where a, ; is the scattering length for interaction between 
molecules in states i and /'. This is isotropic in the case 

a\\ = «22 = an- 
Averaging over the radial and vertical wavefunctions gives 
an effective ID GPE in <f>, 



-fcrot^ + S(0)|O(0)| 2 



<D(0) = 0<D(0, 



(6) 



where 



u(<p) 



4nh 2 Na((f>) 
M 



ff \m\\(Z)\ A pd P dZ (7) 
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and/i = p.+(yz + ^tnjL>z+{v p + ^jtuL> p -b rot /4. For a flat ring, 

this has analytical solutions 0(0) = (27r) -1 ' 2 exp(z'm0), with 
fi = b mt m 2 + u/(2n). Applying antiperiodic boundary condi- 
tions, we obtain the same solutions as for the single-particle 
case, namely states of half-integer quantized angular momen- 
tum. 

In the presence of a small anisotropy, there are two classes 
of solution satisfying antiperiodic boundary conditions: flow- 
ing solutions, <&%((p) ~ (27r)~ 1 ^ 2 exp(±im<f>), and static so- 
lutions, such as » 7r~'/ 2 cos(m0) for small interac- 
tions, both with half-integer m. For a trap with a residual 
anisotropy Vtrap($0 = —V\ cos (/> — Vz cos 2<p and an angle- 
dependent interaction strength u((/>) - uo + ui cos <p + U2 cos 2<f>, 
we obtain approximate chemical potentials corresponding to 
these two Ansatzes of p.f^ 2 » b rot /4 + uo/(2n) and jiy 2 « 
^rot/4 — Vi/2 + 3u()/(4t:) - ui/(2n) - Uzli^Ti). In this approx- 
imation, a flowing state with m = +1/2 is the ground state if 
2m i + U2/2 + 2nV\ < M(>. For any given a(<p) and condensate 
number, we can apply an offset of the optical potential suffi- 
cient to compensate for the anisotropy of the interaction term 
and stabilize the flowing state. 

The permitted velocities of the persistent flow can assume 
only half-integer values compared to the quantized units of 
circulation possible in a more conventional single-species 
atomic BEC in a comparable toroidal geometry [11]. The 
persistent flow should be observable by releasing the trapped 
particles and employing a time-of-flight technique Hill . 

A BEC such as described here is stable only if «(</>) re- 
mains positive all around the ring. From Eq. (Q, this re- 
quires that an and «22 are both positive and that 2ai2 > 
-(an + 022)- In addition, a condensate of polar molecules can 



undergo dipolar collapse [26] if the dipole length a& exceeds 
the scattering length for the short-range interactions, where 
«d = \d 2 \M / (127T€oh 2 ) and d is the effective dipole moment of 
the molecule in the field. Since the molecular wavefunction 
is given by Eq. [4] and there is a direct dipole moment matrix 
element (1|//|2) between the near-degenerate states if/i and 0% 
d = (l|/i|2) sin0. For the two states of KRb considered here, 
{\\H\2) ~ 1(T 32 Cm, which gives a™ x ~5x 1(T 12 m. This is 
substantially smaller than typical scattering lengths so dipolar 
collapse is unlikely. 

The present Letter has described a novel form of conical 
intersection that can occur as function of 3 -dimensional labo- 
ratory coordinates, instead of internal molecular coordinates. 
Using such conical intersections, it may be possible to cre- 
ate novel superfluid states with stable persistent flow charac- 
terized by half-integer, rather than integer quantized angular 
momentum. Although we have considered the effect for po- 
lar molecules, a similar effect might be produced for well- 
separated levels, perhaps even in atoms, using a laser field to 
bring the levels into near-degeneracy and an inhomogeneous 
magnetic field to provide a crossing. Conical intersections 
would appear at points of zero laser amplitude, for example at 
the nodes in an optical lattice. 

The effect proposed here can in principle be observed in 
any system where two states of any different symmetry (not 
just parity) can be tuned into degeneracy with one external 
influence, and then split apart again with another influence 
that breaks the symmetry. 
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